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We present a numerical finite size scaling study of the localization length in long cylinders near 
the integer quantum Hall transition (IQHT) employing the Chalker-Coddington network model. 
Corrections to scaling that decay slowly with increasing system size make this analysis a very 
challenging numerical problem. In this work we develop a novel method of stability analysis that 
allows for a better estimate of error bars. Applying the new method we find consistent results 
when keeping second (or higher) order terms of the leading irrelevant scaling field. The knowledge 
of the associated (negative) irrelevant exponent y is crucial for a precise determination of other 
critical exponents, including multifractal spectra of wave functions. We estimate \y\ > 0.4, which 
is considerably larger than most recently reported values. Within this approach we obtain the 
localization length exponent 2.62 ± 0.06 confirming recent results. Our stability analysis has broad 
applicability to other observables at IQHT, as well as other critical points where corrections to 
scaling are present. 
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Despite a long history, scaling properties of the integer 
quantum Hall transition (IQHT) still pose a consider- 
able challenge [Q. In addition to the critical exponent 
V of the localization length near the transition, recent 
studies (2]-|| analyze the multifractal spectrum A, that 
describes the scaling of moments of the local density of 
states (LDOS) with the system size: {g{r)'^) ^ L^^". 
This spectrum can, in principle, be measured in STM 
experiments. Though promising experiments have been 
undertaken ^, the presently attainable energy reso- 
lution appears to be insufficient to allow a sufficiently 
accurate measurement of Ag. Hence, one relies on nu- 
merical simulations of critical wave functions statistics 
which relate to A, via ((L'^|i/'(r)P)*) - L-^" §. 

On the other hand, one can hope to calculate the spec- 
trum Aq analytically. In the past decade, several pro- 
posals have been made for the quantum field theory un- 
derlying the IQHT critical point |0-|l|] . However, most 
of them rely on assumptions that can neither be taken 
for granted nor easily checked against experiments. The 
identification of the correct critical theory remains an 
outstanding open problem, and predictions of the pro- 
posed models have to be taken with a grain of salt. For 
example, Wcss-Zumino-type theories |f7|-p^ are expected 
to exhibit a strictly parabolic multifractality spectrum 
Aq oc q{l — q). However, the existing numerical evidence 
contradicts this prediction because Aq shows a signifi- 
cant, albeit small, quartic component ^ {q— 1/2)^ ||. 

Generally, numerical studies are performed for finite 
systems. How well the true asymptotic scaling regime 
can actually be probed, often depends crucially on a care- 
ful analysis of subleading corrections to scaling near a 



critical point. They tend to mask the true long-distance 
asymptotics and usually constitute the main difficulty for 
analysis of high-precision numerical simulations. 

Two different kinds of subleading corrections should 
be distinguished. They can be easily understood in the 
framework of renormalization group (RG), and we adopt 
the corresponding terminology in the following discus- 
sion. One kind of subleading behavior is due to the 
fact that expectation values (...) are calculated using a 
Hamiltonian (or an action functional) whose parameters 
(coupling constants) lie on a critical surface, but which 
does not coincide with the fixed point Hamiltonian. De- 
viations from the fixed point along the critical surface are 
parametrized by irrelevant scaling fields, whose decay un- 
der RG flow is described by irrelevant exponents yi < 0. 
The leading irrelevant field decays with the exponent y 
with the smallest absolute value. Its contributions are 
present in subleading corrections to scaling of any corre- 
lation function in a finite system at criticality. 

The second kind of subleading corrections is associated 
with the specific correlation functions that one studies. 
The point is that a particular physical observable of in- 
terest may not exhibit pure scaling behavior even at the 
RG fixed point. Only certain carefully chosen observ- 
ables correspond to pure scaling operators in the fixed 
point theory. As has recently been emphasized, moments 
of the LDOS and critical wave functions at an Anderson 
transition are examples of such pure scaling operators 
]l3[ . Thus, for moments of critical wave functions we ex- 
pect corrections to scaling to come only from irrelevant 
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fields, and primarily from the leading irrelevant field: 
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Even though the correction terms (i^, . . .) do not in- 
fluence scaling to leading order (hence the term "irrel- 
evant"), they are still important to be studied because 
they ultimately determine the size of the critical param- 
eter window (the true scaling regime) iQ. As opposed 
to moments of LDOS, moments of the so-called point 
contact conductances are not pure scaling operators. We 
will return to this issue in subsequent papers. 

For many Anderson transitions corrections to scaling 
decay sufficiently fast, so that the true asymptotic be- 
havior is reliably addressed by a simple scaling analysis 
of numerical simulations [l6|-^. The situation ap- 
pears to be much less favorable for the IQHT. Over time 
it became clear that the Chalker-Coddington (CC) net- 
work model specifically designed for numerical 
analysis, exhibits significant corrections to scaling which 
are very difficult to take into account systematically Q]. 

The problem manifests itself in the recently reported 
low value \y\ « 0.17 |Q, which is smaller than earlier 
estimates by more than a factor of two Q . Recent de- 
velopments climaxed when Amado et al. |^] reported 
the presence of logarithmic corrections near the critical 
point in the CC model, implying, in a sense, y = 0. 
Concomitant with this was an estimate for the ratio 
Tc of the quasi-one dimensional localization length and 
the system width at the IQHT. Conformal invariance, 
which is expected to hold at the IQHT, predicts the 
relation Fc/tt = ao — 2, where = dlS.q/ dq\^^^ ||25[ - 
The values - 2 ^ 0.2596 and 0.2617 obtained 
in Refs. [|[ ^ are drastically different from the estimate 
Fc/tt = 0.223(0.219,0.228] reported in Ref. H. Taken 
at face value, the discrepancy would indicate a break- 
down of conformal invariance at the IQHT fixed point. 
In addition, if the value \y\ were indeed 0.17 or smaller, 
the critical window for the IQHT would hardly be ac- 
cessible with the currently attainable system sizes. As 
a consequence, reliable estimates for other critical expo- 
nents, like Aq, could not be obtained. 

Motivated by this situation we revisit the finite size 
scaling corrections to the localization length in quasi-one 
dimensional geometry (QID) near the IQHT. In addi- 
tion to including the leading irrelevant scaling field in 
our analysis, we develop and employ a new method of 
"stability map" in the parameter space. The fitting pro- 
cedure adjusts two critical exponents (the leading irrel- 
evant exponent y and the exponent of the localization 
length v) and four or more coefficients in the expansion 
of the scaling function. Despite many fitting parameters, 
reliability of our results is demonstrated by the stabil- 
ity map analysis. As a result, we arrive at the following 
conclusions: the estimate of the critical exponent for the 
localization length v = 2.62 ± 0.06 obtained by previous 




FIG. 1. (Color online) (a) rjv/(a;) as a function of for the 
circumference (number of links in the transverse direction) 
M = 16,20,24,32,36,40,48,64,80,96,128,160,192,256,384 
(from the bottom to the top at ~ 0.005). The relative 
standard errors (one sigma) of our data are 0.005% for M < 
64, 0.02% for 80 < M < 192 and 0.05% for M = 256, 384. 
The solid lines are obtained from fitting the first two terms in 
Eq. (^ to the data for each M. Inset: Enlargement of rM{x) 
near x — 0. (b) Same as (a) but in a wider range of x^. 



authors ||15|, ^ g^, |27|-|29|] is confirmed. Furthermore, we 
find \y\ > 0.4 and F^/tt = 0.257 ±0.002. The latter value 
is consistent with predictions coming from conformal in- 
variance, as well as earlier estimates 

To calculate the localization length £,m{x) in QID we 
employ the isotropic version of the CC network model 
pl| , p2| on cylinders (meaning, with periodic boundary 
conditions in one direction) of length L and circumfer- 
ence M. £,m{x) is calculated from the standard trans- 
fer matrix method |2^, The relevant scaling field x 
parametrizes the non-random part of the transfer matrix 



at a node as 
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= 1 - t'^. 

To extract critical exponents we fit numerical data for 
the finite size ratio Fa/(x) = M/£^m{x) by the following 
expression for the scaling function: 



F(x, M) 
7(M) 
7'(M) 



7(M) + x^j'{M) + Oix'^M'^'"), 



(2) 
(3) 
(4) 



(Notice that we distinguish the data Fjv/(a:;) from the scal- 
ing function T{x,M) by using different notation.) The 
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FIG. 2. (Color online) Original data Tm{x = 0) (black x, 
left Eixis) and the offset 7m (red o, left axis) and the slope 'j'^j 
(blue •, right axis) obtained ftom fitting the first two terms 
in Eq. (||) to Fa/ (a:) for x < 0.05. 

analytical form of Eqs. (j^-^ is standard for a scaling 
function near a critical fixed point combined with a fact 
that r{x,M) is an even function of x at the IQHT [ p3| . 
It involves an expansion in even powers of x and integer 
powers of the leading irrelevant field that itself is assumed 
to scale as with the systems width. We keep only the 
leading order in x^M^^'^ to reduce the number of fitting 
parameters. Moreover, while we will keep the terms up 
to for 7(M), we will ignore terms starting from M^^ 
for 7'(M). The reason is that within our range of system 
sizes keeping these correction terms changes the exponent 
ly at the level of less than one percent. 

We show numerical data Fa./ (cc) in Fig. |^. The solid 
straight lines in the figure show the results of the least 
square fit of the first two terms in Eq. (||) to the data, 
performed separately for each M in the range x < 0.05. 
We will denote the offset and the slope of the straight 
lines obtained in this way by jm and 7^^. We see from the 
plots that corrections coming from cc^ (and higher order 
terms) are very small even within the range of x'^ shown 
in Fig. 0(b). This justifies keeping only terms up to order 
x^ in the expansion (||). The effects of corrections to 
scaling are most pronounced near the critical point (Fig. 
0, inset) where the solid lines intersect at nonzero values 
of X instead of meeting at a single point Fc at a; = 0. 

In Fig. H we plot the offset and the slope 7^/ ob- 
tained from fitting by Eq. (|^) , together with the original 
data for Tm^x — 0). We see that the all these quantities 
significantly depend on the system size M . The variation 
of Fa/(0) over the available range of M is ^ 4%. The very 
slow (approximately logarithmic) dependence on M seen 
in Fig. ^ constitutes the notorious difficulty for numeri- 
cal studies on the IQHT mentioned in the introduction. 
At this point we can proceed in two different ways. 

Method I (the two-step optimization method). In this 
method we take a subset of Fm(0) and 7^ presented in 
Fig. I and use these as the input data set for fitting 
the functions 7(M) and -/'{M) from Eqs. (|, |). These 
functions depend on a number of fitting parameters v, 
y, Fc, F', ai, a'l, etc., that we optimize by finding a 



minimum to the "cost" function 
^ (Fm(0)- 7(M))^ 



(0) 



M 



'M 



where aM{^),cr'M ^'^s the (absolute) statistical errors of 
the input data Fm(0), 7^^ Q. 

Method II (the global optimization method). In this 
method we take a subset of all numerical data Fjv/(a;) as 
the input for a global fitting by the function F(a;, M) from 
Eqs. (||-^). Optimal parameters of F(x, M) are found by 
minimizing the cost function 
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(FM(x)-F(a;,M))= 



'M 



(x) 



(6) 



where (Jm{x) are the (absolute) standard errors of FmC^;). 

In both methods we can vary the input sets by varying 
the range Mmin-Mmax of the system sizes and the total 
number N of data points in the set. 

We have employed both methods I and II with more 
than 1000 fitting trials for each input data set, starting 
with different initial fitting parameters chosen randomly 
p2[ . To visualize the results, we show a "stability map" 
in Fig. ^. It displays minima of Xi n functions of all 
fitting parameters. The map allows for a better estimate 
of error bars and helps understanding how fitting param- 
eters affect the cost functions. 

For method II the optimization routine yields a unique 
minimum of Xn for each data set. The minima for 
sets with different ranges of M and different numbers 
of points N are tightly clustered together if we take into 
account the terms up to a2M^^ in the function 7(M), but 
not when we keep only the first correction term aiM^, 
see lower panels in Fig. ||. 

In stark contrast, the cost function Xj exhibits multiple 
local minima where method I gets stuck. On the upper 
panels in Fig. |^ we show the local minima with low values 
of Xi obtained by method I with terms up to a2M'^y m 
7(M). From the panel (Tb) one infers that the local 
minima exist in a wide range Q < \y\ < 0.7. However, 
only for values of \y\ above \y\ k, 0.4 the cost function Xi 
approaches its global minimum, which is emphasized by 
colored symbols with error bars. At this global minimum 
the results of methods I and II agree. 

From the fitted data presented in Fig. |^ we extract 



several observations 1 33 



(i) If only the aiM^ term is included in 7(M), and 
Mjnin < 40, we find \y\ < 0.2. This is consistent with a 
previously reported value [Q. The small apparent ex- 
ponent y reflects the slow decay of 7m with 1/M seen 
in Fig. 1^. Remarkably, when the next correction term 
is kept in 7(Af), the optimal value for \y\ increases; we 
obtain 0.4 < \y\ < 0.6 |35j. This estimate remains un- 
changed when we add the term osM'^^'lQ. Larger values 
of \y\ obtained when ai^2-terms and higher ones are in- 
cluded do not contradict the slow decay of ff these 
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with a2 in y(Af) 
16-384(30) 

16-256 28) I 

20-384 28) ^-m 
20-256(26) I — E 
24-384(26) i — * 
24-256(24) i ^ 
32-384(24) i — » 
32-256 22) ^ 
36-384 22) i — • 
36-256(20) I — 



with 02 in y(M) 
16-384(176) 
16-256 168 
20-384 161) ^-m- 
20-256 153 ^a- 
24-384(146) i — *— 
24-256(138) i — ^ 
32-384(131) I — T— 
32-256 123 
36-384(116) I — 
36-256(108) ^9— 
without 02 in y{M) 
32-384(131) I — 
32-256(123) i — X— 
36-384(116) I — 
36-256 108 ^s— 
40-384(101) I — ©— 
40-256(93 ^e— 
48-384(86) i — 9— 
48-256(78) i — e— 
64-384(71) I — G— 
64-256(63) i e 



FIG. 3. (Color online) Stability map for the exponents u, \y\ and the coefficients Fc/tt, ai, 02, F', a'l appearing in Eqs. 
(i) obtained by methods I and II for < x < 0.05. In a given panel each symbol shows a local minimum of the cost 
function Xi.ii divided by the number A'' of data points taken into account. Different symbols represent different input data sets 
obtained by varying the range of system sizes M and the number of data points A*', and listed in the legend boxes at the right 
as Mmin-M^s,^{N). Upper panel: In the case of Xi many local minima (shown by gray symbols) have been found. The global 
minima of Xi for each input data set are shown by colored symbols with error bars |^^. Lower panel: For xii a single minimum 
(colored symbols) was found for each data set, when the term a2M^^ in Eq. (^ was included in the analysis. We also show 
results of minimization (black symbols) for the case when a single correction term ailVP in 7(M) was taken into account for 
input data with Mmin < 64. On the panel (Il-e) the black symbols are shown on the vertical line 02 = 0. We note that since 
/N tends to be small for small A'', a subtle difference of x^ between different input data sets is meaningless. 



terms partially cancel each other for small system sizes. 
Indeed, we see from Fig. || that the coefficients ai and 
a2 are similar in magnitude but opposite in sign. 

(ii) The estimate of Fc/tt = ag ~ 2 is sensitive to the 
number of correction terms kept in 7(M). If only the ai 
term is kept, and A/min < 40, our value of Fc is again 
consistent with the result ao - 2 « 0.248[0.244, 0.251] 
found earlier |^3|. However, adding higher order terms 
gives a significantly larger value ao — 2 « 0.257 ± 0.002. 
This latter estimate is broadly compatible with results 
based on the wave function statistics [H-Q. 

(iii) We confirm the estimate for the localization length 
exponent 2.62 ± 0.0 6 a lready obtained previously by 
several authors p4|, p7|-p9|. 

Given that finite size effects are quite significant near 
the IQHT, and that they have been treated very differ- 
ently by previous authors, it is, perhaps, remarkable that 
similar estimates for v were recently found. Using the 
same model as in this work, Ref. kept only a sin- 
gle My term, and Ref. |2^ employed powers of l/ln7\/. 
Different boundary conditions or models were employed 
in Refs. |2^, |2^, where only a single term was used. 
The authors of Ref. |2^ found corrections to scaling to be 
insignificant. We believe that this relative insensitivity 
of the estimates for v to subleading corrections is related 



to the facts that (a) v is large compared to \y\ and (b) 
the ranges of system sizes M used are narrow, with the 
ratio Afniax/-Afmin hardly exceeding ten. 

By the same reasoning, it is not clear why earlier 
estimates of v obtained for the random-Landau-matrix 
model [v = 2.35 ± 0.03 13) or other models [v = 
2.33 ±0.09 m, ^), differ from more recent values by up 
to 10%. While a violation of universality at the IQHT 
would be a logically possible explanation, the present 
numerical evidence is not strong enough to draw such 
a drastic conclusion. We expect that the apparent dis- 
crepancy in estimates of scaling exponents obtained from 
different models of the IQHT will disappear upon rein- 
vestigating finite size effects more carefully. 

Conclusions. Using the CC network model on long 
cylinders, we have numerically analyzed corrections to 
scaling near the IQHT. Our data is consistently inter- 
preted using the standard form of corrections to scaling 
if more than one correction term is included in the scaling 
functions. The interpretation is facilitated by the "sta- 
bility map" analysis. Our results satisfy predictions of 
conformal invariance. 

We thank S. Bera, A. Furusaki, A. D. Mirlin, T. Oht- 
suki, A. Sedrakyan, K. Slevin and A. Tsvelik for numer- 
ous discussions of corrections to scaling. We thank K. 



5 



Slcvin and T. Ohtsuki for sharing Rcf. [|T^ prior to pub- 
lication, and A. W. W. Ludwig for helpful comments. We 
also thank I. Kondov (SCC) for support and the Juelich 
Supercomputer Center (project HKA12) for allocation 
of computing time on JUROPA. Numerical simulations 
were performed in part using the PADS resource (NSF 
Grant No. OCI-0821678) at the Computational Institute 
at the University of Chicago. H. O. is supported by a 
Grant-in-Aid for Research Abroad from JSPS. I. A. G. 
was supported by NSF Grants No. DMR-1105509 and 
No. DMR-0820054. 



* Present address: Department of Applied Physics, 
Hokkaido University, Sapporo 060-8628, Japan 
[1] F. Evers and A. D. Mirlin, Rev. Mod. Phys. 80, 1355 
(2008). 

[2] F. Evers, A. Mildenberger, and A. D. Mirlin, Phys. Rev. 

B 64, 241303(R) (2001). 
[3] H. Obuse, A. R. Subramaniam, A. Furusaki, I. A. 

Gruzberg, and A. W. W. Ludwig, Phys. Rev. Lett. 101, 

116802 (2008). 

[4] F. Evers, A. Mildenberger, and A. D. Mirlin, Phys. Rev. 

Lett. 101, 116803 (2008). 
[5] M. Morgenstern, J. Klijn, C. Meyer, and R. Wiesendan- 

ger, Phys. Rev. Lett. 90, 056804 (2003). 
[6] K. Hashimoto, C. Sohrmann, J. Wiebe, T. Inaoka, F. 

Meier, Y. Hirayama, R. A. Romer, R. Wiesendanger, and 

M. Morgenstern, Phys. Rev. Lett. 101, 256802 (2008). 
[7] M. R. Zirnbauer, arXiv:hep-th/9905054. 
[8] M. J. Bhaseen, I. I. Kogan, O. A. Soloviev, N. Taniguchi, 

and A. M. Tsvelik, Nucl. Phys. B580, 688 (2000). 
[9] A. M. Tsvelik, Phys. Rev. B 75, 184201 (2007). 
[10] A. LeClair, arXiv:0710.3778vl. 

[11] C. A. Liitken, and G. G. Ross, Phys. Lett. A 356, 382 
(2006); Phys. Lett. B 653, 363 (2007); Nucl. Phys. B850, 
321 (2011). 

[12] E. Bettelheim, I. A. Gruzberg, and A. W. W. Ludwig, 
Phys. Rev. B 86, 165324 (2012). 

[13] I. A. Gruzberg, A. W. W. Ludwig, A. D. Mirlin, and M. 
R. Zirnbauer, Phys. Rev. Lett. 107, 086403 (2011). 

[14] Irrelevant corrections to scaling play an important role 
in the work by Liitken and Ross [^]. These authors 
postulate an anti-holomorphic symmetry property of the 
scaling function I3^=l3xy+il3xx in the coupling constant 
o'=axy+icrxx, and use it to derive the relation v = 
10.42 . . . /c where c is an unspecified integer [ ]ll] , [l^ ]. The 
choice c = 4 would lead to values consistent with numeri- 
cal estimates. A necessary condition for anti-holomorphic 
scaling is = 1. This condition is not inconsistent with 
our data in Fig. I (IIa,IIb), since the colored symbols 
correspond to values of \y\u between 1.0 and 1.5 with 
pronounced clustering in the range 1.0-1.3. For this pos- 
sibility, see also Ref. hs. 

[15] K. Slevin and T. Ohtsuki, Int. J. Mod. Phys. Conf. Ser. 
11, 60 (2012). 

[16] K. Slevin and T. Ohtsuki, Phys. Rev. Lett. 82, 382 
(1999). 

[17] V. Kagalovsky, B. Horovitz, Y. Avishai, and J. T. 
Chalker, Phys. Rev. Lett. 82, 3516 (1999). 



[18] Y. Asada, K. Slevin, and T. Ohtsuki, Phys. Rev. Lett. 

89, 256601 (2002); Phys. Rev. B 70, 035115 (2004). 
[19] A. Rodriguez, L. J. Vasquez, K. Slevin, and R. A. Romer, 

Phys. Rev. Lett. 105, 046403 (2010); Phys. Rev. B 84, 

134209 (2011). 

[20] R. Bondesan, I. A. Gruzberg, J. L. Jacobsen, H. Obuse, 
and H. Saleur, Phys. Rev. Lett. 108, 126801 (2012). 

[21] J. T. Chalker and P. D. Coddington, J. Phys. C, 21, 2665 
(1988). 

[22] B. Kramer, T. Ohtsuki, and S. Kettemann, Phys. Rep. 

417, 211 (2005). 
[23] K. Slevin and T. Ohtsuki, Phys. Rev. B 80, 041304(R) 

(2009). 

[24] M. Amado, A. V. Malyshev, A. Sedrakyan, and F. 

Dommguez-Adame, Phys. Rev. Lett. 107, 066402 (2011). 
[25] M. Janssen, Int. J. Phys. B 8, 943 (1994); Phys. Rep. 

295, 1 (1998). 

[26] A. Dohmen, P. Freche, and M. Janssen, Phys. Rev. Lett. 

76, 4207 (1996). 
[27] H. Obuse, A. R. Subramaniam, A. Furusaki, I. A. 

Gruzberg, and A. W. W. Ludwig, Phys. Rev. B 82, 

035309 (2010). 

[28] J. P. Dahlhaus, J. M. Edge, J. Tworzydlo, and C. W. J. 
Beenakker, Phys. Rev. B 84, 115133 (2011). 

[29] I. C. Fulga, F. Hassler, A. R. Akhmerov, and C. W. J. 
Beenakker, Phys. Rev. B 84, 245447 (2011). 

[30] A. MacKinnon and B. Kramer, Z. Phys. B: Condens. 
Matter 53, 1 (1983). 

[31] Alternatively, we could have used the fitted data 7i\/ in- 
stead of rA/(0) for finding the optimal parameters of the 
functions 7(Af) and 7'(Af) (by modifying the cost func- 
tion appropriately). Since both data sets, rj\/(0) and 7a/, 
hardly deviate from each other (see Fig. both proce- 
dures are equivalent in the present case. 

[32] We have given trivial constraints to the fitting: y is neg- 
ative, contributions of correction terms become smaller 
than the main terms, Fc and F', at the maximal M, and 
the goodness of fit is larger than 0.01. While we mainly 
discuss the case when x from the interval < a; < 0.05 
(F < 1) was used in the fitting, we have checked that our 
fits are stable against varying the range of 2;-values that 
we use js^ ]. 

[33] See Supplemental Material for details of numerical simu- 
lations. There we also provide a comparison of our resutls 
with previously published ones. 

[34] The estimation of error bars for fitting parameters is 
based on propagation of uncertainty. Therefore, these er- 
ror bars represent volatility against the statistical fluctu- 
ation of raw data rA/(a;) within its error bar aM{x). 

[35] Moreover, in method I, at the local minima 

of X? 

is almost zero, and the corresponding values of \y\ are 
located in \y\ < 0.2, and Fc/tt < 0.25, as shown in Fig. 
I (Lb). This is consistent with the results of method II, 
when only the ai term is kept in 7(M). 
[36] We estimate u based on data obtained by method II in- 
cluding second order correction terms in MM). As can 
be seen from Fig. ^ (Ha) and also in Ref. [^ (Tables III 
and IV) there is a small residual scattering of the fitting 
parameters v with variation of the input data. Our es- 
timate V = 2.62 ± 0.06 exhibits error bars that describe 
the full width of the window of i^-values that we have 
observed. 

[37] B. Huckestein, Rev. Mod. Phys. 67, 357 (1995) and ref- 
erences therein. 
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[38] N. Sandler, H. R. Maei, and J. Kondev, Phys. Rev. B 70, 

045309 (2004). 
[39] See also references in Ref. ^ 
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Supplemental Material for "Finite Size 
Effects and Irrelevant Corrections to Scal- 
ing near the Integer Quantum Hall Tran- 
sition" 

In this supplemental material, we present all fitting pa- 
rameters for scaling functions with up to three correction 
terms, for various input data sets. We also show the 
stability maps for these parameters, and provide a table 
summarizing and comparing our results for critical expo- 
nents with those of other publications. 



SCALING ANALYSIS 

We calculate the localization length for the 

Chalker-Coddington network model in the quasi-one di- 
mensional cylinder geometry close to the integer quan- 
tum Hall transition (at x = 0). We use the following 
values for the circumference of the cylinder M = 16, 20, 
24, 32, 36, 40, 48, 64, 80, 96, 128, 160, 192, 256, 384. To 
investigate the critical behavior, we employ the following 
scaling function for Tn[{x) = M/^m{x), 

T{x,M)=^{M)+x^i{M), (1) 
7(A/) - r, (1 + aiAP + a2M^y + ogM^y) , (2) 
7'(A/) = r'A/2/>^ (1 + a'^Np) , (3) 

where the fitting parameters v and y represent the critical 
exponent of the localization length and the leading irrel- 
evant exponent, respectively. The other fitting parame- 
ters are the coefficients Fc, ai, 02, 03, F', a'l of the scaling 
functions. 

As is explained in the main paper, we consider two 
fitting methods. In method I we use the data set of 
Fa/(.t = 0), and 7^^ obtained from fitting Eq. to 
the data for Tf,[{x). In method II we use the original 
data T m{x). Furthermore, in order to check that our re- 
sults remain unchanged upon variations in the x range 
used in the fitting, we consider two ranges: < a; < 
with Xu — 0.05 and x^ ~ 0.07. Thereby, we have four 
kinds of input data sets total. 

For each method, we seek the global minimum of the 
corresponding cost function ^ (see Eqs. (5-6) in the 
main paper). We examine the stability of fitting with 
respect to keeping different number of correction terms 
in in Eq. (|^), from one to three. We also vary the min- 
imum and maximum circumference, Afniin and .A /max, of 
the input data sets to check stability of our fits. 

Tables | - IV shows results of fitting corresponding to 
the global minimum of for the four kinds of input data, 
i.e., method I with x^ = 0.05, method I with x„ = 0.07, 
method II with x„ — 0.05, and method II with x„ = 0.07. 
In addition, we show the stability maps for these input 
data in Fig. |[ 

Some remarks are in order. 



(i) Even if the input data for Xu = 0.07 are used, fit- 
ting parameters for the global minimum are consistent 
with those obtained from the input data for x^ = 0.05. 
Moreover, distributions of the local minima obtained by 
method I as shown in Fig. |l] are also similar. For exam- 
ple, the local minima for 02 locate near 02 = 0, and most 
of those for \y\ are located below \y\ = 0.2, etc. This 
demonstrates stability of fitting procedure with respect 
to changes in x^. 

(ii) In addition to the first two correction terms, we 
have examined the effect of the third correction term 
a^M^y in the scaling function 7(M). Since the higher 
order terms can contribute appreciably to Fm (x) only for 
small M, we fix Afmin to a small value, 16. (For larger 
values of A/min results of the fitting become unreliable, 
which manifests itself in error bars of fitting parameters 
exceeding their mean values.) Keeping the term a^M^y 
in the sealing function results in changes of all fitting pa- 
rameters by only a few percent. In contrast, when we go 
from one correction term in the scaling function to two 
(up to 02 A/^^), some fitting parameters change drasti- 
cally (for example, \y\ increases by a factor of 3, the co- 
efficient ai grows by a factor of 2, etc.). We also observe 
that the signs of oi and 02 are opposite when including 
the second correction term, and this feature remains un- 
changed when adding the third correction term to the 
analysis. As we have explained in the main paper, the 
opposite signs lead to the estimate of |y| to grow due to 
partial cancellations. These results imply that keeping 
two correction terms in 7 (A/) is sufficient for our pur- 
poses. 



COMPARISON WITH PREVIOUS WORKS 

For readers' convenience, we compare our results for 
V, y, and Fc/tt (= ao — 2) with those of other authors in 
Table 0. 
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FIG. 1. Stability maps for the exponents v, \y\ and the coefficients Fc/tt, oi, 02, F', a'l. From top to bottom: method I with 
Xn = 0.05 [(Fa)-(Fg)], method I with Xy. = 0.07 [(Fh)-(I-n)], method II with Xu = 0.05 [(Il-a)-(II-g)], and method II with 
a;„ =0.07 [(Il-h)-(II-n)]. 



TABLE I. List of fitting parameters corresponding to the global minimum Xmin of the cost function, obtained by method I for < a; < 0.05. i\/min and Mmax represent 
the minimum and maximum circumferences, respectively, used in the analysis. and Q indicate the total number of input data points and goodness of fit, respectively. 



A/min A/max A'' 


Xmin/iV Q 


u \y\ Fc/tt = ao - d 


ai a') a-j 


1 ai 


24 384 26 
32 384 24 
36 384 22 
40 384 20 

48 '^S4 1 8 

64 384 16 


1.008 0.159 
0.935 0.213 
0.828 0.311 
0.680 0.480 

n fi4^ fl 477 

0.497 0.634 


2.712 ± 0.055 0.136 ± 0.012 0.24349 ± 0.00203 
2.669 ± 0.059 0.162 ± 0.020 0.24699 ± 0.00231 
2.622 ± 0.062 0.181 ± 0.024 0.24886 ± 0.00225 
2.567 ± 0.067 0.200 ± 0.028 0.25035 ± 0.00210 

9 '^dR + n flfiS n 94'^ + fl 04^ fl 9'i981 + fl 0091 fl 

2.651 ± 0.086 0.255 ± 0.064 0.25333 ± 0.00274 


0.1639 ±0.0078 
0.1519 ±0.0073 
0.1470 ±0.0058 
0.1442 ±0.0039 

n 144^1 + flnn9^ — — 
0.1457 ±0.0067 


0.768 ±0.139 -0.442 ±0.203 
0.657 ±0.120 -0.329 ±0.214 
0.568 ±0.104 -0.190 ±0.227 
0.481 ± 0.094 0.009 ± 0.255 
fl 4^17 + n fls^ n flss + n 97fl 

0.599 ±0.133 -0.308 ±0.325 


24 256 24 
32 256 22 
36 256 20 
40 256 18 
48 256 16 
64 256 14 


1.017 0.142 
0.993 0.148 
0.908 0.199 
0.745 0.340 
0.653 0.402 
0.427 0.650 


2.749 ± 0.070 0.132 ± 0.013 0.24270 ± 0.00234 
2.701 ± 0.080 0.158 ± 0.023 0.24649 ± 0.00277 
2.635 ± 0.087 0.181 ± 0.029 0.24883 ± 0.00265 
2.544 ± 0.098 0.205 ± 0.034 0.25076 ± 0.00240 
2.500 ± 0.102 0.273 ± 0.054 0.25415 ± 0.00212 
2.718 ±0.142 0.333 ±0.093 0.25596 ± 0.00235 


0.1670 ±0.0093 
0.1536 ±0.0091 
0.1471 ±0.0069 
0.1434 ±0.0041 

0.1463 ± 0.0062 - — 
0.1588 ±0.0237 


0.851 ±0.195 -0.512 ±0.254 
0.717 ±0.175 -0.403 ±0.282 
0.589 ±0.148 -0.228 ±0.308 
0.450 ± 0.129 0.098 ± 0.372 
0.404 ±0.110 0.302 ±0.415 
0.671 ±0.206 -0.557 ±0.531 


24 192 22 
32 192 20 
36 192 18 
40 192 16 
48 192 14 
64 192 12 


1.020 0.130 
1.027 0.114 
0.957 0.141 
0.785 0.249 
0.697 0.283 
0.492 0.434 


2.755 ±0.076 0.127 ±0.013 0.24185 ± 0.00258 
2.705 ± 0.090 0.150 ± 0.024 0.24546 ± 0.00322 
2.630 ± 0.099 0.171 ± 0.030 0.24791 ± 0.00312 
2.516 ± 0.117 0.196 ± 0.036 0.25006 ± 0.00281 
2.451 ± 0.126 0.264 ± 0.059 0.25378 ± 0.00247 
2.739 ± 0.184 0.325 ± 0.104 0.25574 ± 0.00279 


0.1705 ±0.0105 
0.1572 ±0.0112 
0.1497 ±0.0090 

0.1448 ±0.0057 - - 
0.1456 ±0.0053 
0.1568 ±0.0247 


0.873 ±0.221 -0.525 ±0.278 
0.734 ±0.205 -0.416 ±0.320 
0.584 ±0.173 -0.216 ±0.356 
0.413 ±0.147 0.217 ±0.453 
0.350 ±0.127 0.548 ±0.524 
0.702 ± 0.277 -0.609 ± 0.662 


16 384 30 
20 384 28 

94 "^84 9fi 
32 384 24 
36 384 22 


0.875 0.289 
0.813 0.356 

n 7'iQ fl 41 1 

0.755 0.381 
0.724 0.387 


2.655 ± 0.016 0.413 ± 0.027 0.25607 ± 0.00066 
2.640 ± 0.019 0.404 ± 0.040 0.25587 ± 0.00096 
9 fi'Hi + n 09^1 n 'Hfis + fl n^7 fl 9'\'\n7 + n nni 48 

2.611 ±0.025 0.428 ±0.076 0.25623 ± 0.00143 
2.591 ± 0.027 0.450 ± 0.101 0.25656 ± 0.00164 


0.2687 ±0.0116 -0.2843 ±0.0417 
0.2642 ±0.0189 -0.2688 ± 0.0648 

n 944^ + fl n9^8 — n 9fl79 + fl 0784 — 

0.2863 ± 0.0525 -0.3417 ± 0.1904 
0.3050 ± 0.0824 -0.4113 ± 0.3231 


0.579 ±0.018 -0.387 ±0.061 
0.563 ± 0.022 -0.326 ± 0.074 
fl + n fl'Hfl — n 9fi7 + n flQ9 

\J.OOO ZIZ U.UOU KJ.j^UI zn u.uy^ 

0.529 ±0.028 -0.213 ±0.113 
0.508 ±0.028 -0.123 ±0.133 


16 256 28 
20 256 26 
24 256 24 
32 256 22 
36 256 20 


0.810 0.362 
0.783 0.374 
770 359 
0.744 0.358 
0.693 0.384 


2.667 ± 0.019 0.425 ± 0.029 0.25640 ± 0.00068 
2.650 ± 0.023 0.427 ± 0.045 0.25644 ± 0.00098 
2 640 ± 029 394 ± 066 25576 ± 00155 
2.615 ± 0.029 0.496 ± 0.085 0.25743 ± 0.00126 
2.590 ± 0.029 0.569 ± 0.106 0.25826 ± 0.00118 


0.2733 ±0.0129 -0.3029 ± 0.0479 
0.2746 ± 0.0228 -0.3072 ± 0.0825 
0.2558 ± 0.0322 —0.2442 ± 0.1046 — 
0.3331 ±0.0710 -0.5371 ± 0.3224 
0.4147 ±0.1256 -0.9664 ± 0.7638 


0.589 ± 0.021 -0.430 ± 0.069 
0.571 ±0.025 -0.369 ±0.085 
563 ± 035 —0 310 ± 109 
0.531 ±0.030 -0.266 ±0.145 
0.504 ±0.027 -0.154 ±0.186 


16 192 26 
20 192 24 
24 192 22 
32 192 20 
36 192 18 


0.859 0.268 
0.838 0.268 
0.824 0.256 
0.814 0.234 
0.766 0.245 


2.670 ± 0.020 0.422 ± 0.031 0.25631 ± 0.00073 
2.653 ± 0.024 0.421 ± 0.049 0.25628 ± 0.00109 
2.644 ± 0.033 0.380 ± 0.074 0.25538 ± 0.00187 
2.616 ± 0.032 0.485 ± 0.094 0.25726 ± 0.00147 
2.587 ± 0.032 0.564 ± 0.117 0.25820 ± 0.00134 


0.2721 ±0.0131 -0.2976 ± 0.0487 
0.2716 ±0.0235 -0.2959 ± 0.0844 
0.2496 ± 0.0333 -0.2235 ± 0.1059 
0.3254 ± 0.0748 -0.5019 ± 0.3295 
0.4090 ±0.1343 -0.9314 ± 0.8035 


0.593 ± 0.023 -0.436 ± 0.072 
0.574 ± 0.028 -0.374 ± 0.089 
0.568 ±0.040 -0.314 ±0.116 
0.532 ±0.033 -0.262 ±0.153 
0.502 ±0.030 -0.137 ±0.197 


16 384 30 
16 256 28 
16 192 26 


0.874 0.242 
0.801 0.318 
0.854 0.223 


2.651 ± 0.036 0.429 ± 0.154 0.25629 ± 0.00210 
2.653 ±0.031 0.490 ±0.133 0.25722 ± 0.00171 
2.657 ± 0.038 0.475 ± 0.159 0.25701 ± 0.00212 


0.2871 ±0.1659 -0.3668 ± 0.7741 0.0878 ± 0.8625 
0.3437 ±0.1457 -0.6670 ± 0.8674 0.4580 ± 1.2776 
0.3291 ± 0.1695 -0.5825 ± 0.9517 0.3410 ± 1.3033 


0.573 ± 0.050 -0.388 ± 0.068 
0.570 ± 0.038 -0.436 ± 0.070 
0.575 ± 0.049 -0.439 ± 0.074 



TABLE II. List of fitting parameters corresponding 
represent the minimum and maximum circumferences, 
respectively. 



to the global minimum Xmin of the cost function, obtained by method I for < a; < 0.07. A/min and Mmax 
respectively, used in the analysis. A'^ and Q indicate the total number of input data points and goodness of fit. 



iV 



\y\ 



Fc/tt = ao — d 



ai 



"2 



"3 



r' 



24 
32 
36 
40 
48 
64 



384 
384 
384 
384 
384 
384 



1.075 
0.868 
0.855 
0.817 
0.647 
0.722 



0.111 
0.288 
0.278 
0.293 
0.474 
0.316 



2.669 ± 0.027 
2.628 ± 0.028 
2.618 ±0.029 
2.604 ± 0.030 
2.566 ±0.031 
2.570 ± 0.040 



0.137 
0.163 
0.182 
0.200 
0.243 
0.255 



±0.012 
± 0.020 
± 0.024 
± 0.028 
± 0.043 
± 0.064 



0.24354 
0.24702 
0.24891 
0.25039 
0.25281 
0.25333 



± 0.00202 
± 0.00230 
± 0.00223 
± 0.00209 
± 0.00210 
± 0.00275 



0.1637 
0.1518 
0.1468 
0.1441 
0.1445 
0.1457 



± 0.0078 
± 0.0072 
± 0.0057 
± 0.0039 
± 0.0025 
± 0.0067 



0.692 ± 0.066 
0.597 ± 0.056 
0.574 ± 0.052 
0.547 ± 0.049 
0.489 ± 0.041 
0.495 ± 0.053 



-0.369 ±0.105 
-0.249 ±0.107 
-0.218 ±0.108 
-0.170 ±0.113 
-0.034 ±0.121 
-0.051 ±0.161 



24 
32 
36 
40 
48 
64 



256 
256 
256 
256 
256 
256 



0.817 
0.729 
0.707 
0.674 
0.558 
0.446 



0.355 
0.450 
0.439 
0.435 
0.539 
0.620 



2.716 ±0.034 
2.675 ± 0.039 
2.670 ± 0.041 
2.659 ± 0.044 
2.612 ±0.046 
2.672 ± 0.070 



0.132 
0.158 
0.181 
0.205 
0.272 
0.334 



±0.013 
± 0.023 
± 0.029 
± 0.034 
± 0.054 
± 0.093 



0.24271 
0.24647 
0.24882 
0.25076 
0.25410 
0.25599 



± 0.00234 
± 0.00278 
± 0.00265 
± 0.00240 
± 0.00215 
± 0.00233 



0.1670 
0.1537 
0.1471 
0.1434 
0.1462 
0.1591 



± 0.0093 
±0.0091 
± 0.0069 
± 0.0041 
± 0.0060 
± 0.0239 



0.793 ± 0.098 
0.682 ± 0.089 
0.657 ± 0.084 
0.627 ± 0.079 
0.546 ± 0.062 
0.616 ±0.100 



-0.469 ±0.133 
-0.368 ±0.145 
-0.352 ±0.147 
-0.323 ±0.154 
-0.194 ±0.168 
-0.433 ±0.258 



24 
32 
36 
40 
48 
64 



192 
192 
192 
192 
192 
192 



0.788 
0.731 
0.728 
0.713 
0.624 
0.447 



0.364 
0.405 
0.362 
0.327 
0.366 
0.498 



2.725 ± 0.037 
2.684 ± 0.044 
2.680 ± 0.047 
2.669 ± 0.050 
2.617 ±0.053 
2.717 ±0.094 



0.127 
0.150 
0.171 
0.196 
0.262 
0.325 



±0.013 
± 0.024 
± 0.030 
± 0.036 
± 0.059 
±0.104 



0.24185 
0.24543 
0.24789 
0.25004 
0.25370 
0.25574 



± 0.00258 
± 0.00323 
± 0.00313 
± 0.00282 
± 0.00250 
± 0.00279 



0.1705 
0.1573 
0.1497 
0.1449 
0.1454 
0.1569 



±0.0105 
±0.0112 
± 0.0090 
± 0.0058 
±0.0051 
±0.0247 



0.819 ±0.112 
0.705 ±0.106 
0.681 ±0.102 
0.648 ± 0.096 
0.554 ± 0.075 
0.676 ±0.149 



-0.488 ±0.146 
-0.391 ±0.166 
-0.377 ±0.169 
-0.349 ±0.176 
-0.209 ±0.193 
-0.558 ±0.322 



16 
20 
24 
32 
36 



384 
384 
384 
384 
384 



1.352 
0.935 
0.790 
0.674 
0.728 



0.013 
0.200 
0.363 
0.512 
0.381 



2.623 ± 0.008 
2.609 ±0.010 
2.603 ±0.013 
2.585 ±0.012 
2.583 ±0.013 



0.419 
0.410 
0.374 
0.432 
0.459 



± 0.027 
± 0.040 
± 0.056 
± 0.075 
± 0.098 



0.25624 
0.25602 
0.25523 
0.25629 
0.25671 



± 0.00063 
± 0.00092 
± 0.00142 
± 0.00140 
± 0.00154 



0.2716 
0.2672 
0.2473 
0.2886 
0.3129 



±0.0117 
±0.0191 
± 0.0260 
± 0.0526 
± 0.0830 



-0.2951 ±0, 
-0.2790 ± 0, 
-0.2158 ±0, 
-0.3501 ±0, 
-0.4426 ± 0, 



0427 
0663 
0799 
1930 
3374 



0.549 ± 0.010 
0.534 ± 0.012 
0.529 ±0.016 
0.508 ±0.014 
0.505 ±0.014 



-0.311 ±0.029 
-0.250 ±0.035 
-0.206 ± 0.044 
-0.151 ±0.053 
-0.149 ±0.062 



16 
20 
24 
32 
36 



256 
256 
256 
256 
256 



0.909 
0.628 
0.569 
0.493 
0.496 



0.228 
0.635 
0.692 
0.763 
0.701 



2.638 ±0.010 
2.623 ±0.012 
2.619 ±0.016 
2.598 ±0.015 
2.595 ±0.015 



0.429 
0.429 
0.394 
0.492 
0.568 



± 0.029 
± 0.045 
± 0.066 
± 0.086 
±0.106 



0.25648 
0.25647 
0.25576 
0.25737 
0.25825 



± 0.00067 
± 0.00097 
± 0.00155 
± 0.00129 
±0.00119 



0.2748 
0.2753 
0.2558 
0.3301 
0.4128 



±0.0130 
± 0.0228 
±0.0321 
±0.0707 
±0.1252 



-0.3085 ± 0, 
-0.3098 ± 0, 
-0.2443 ± 0, 
-0.5237 ±0, 
-0.9547 ± 0, 



0484 
0829 
1044 
3169 
7565 



0.562 ±0.012 
0.546 ± 0.014 
0.544 ± 0.021 
0.518 ±0.016 
0.514 ±0.015 



-0.360 ± 0.033 
-0.302 ± 0.040 
-0.267 ±0.052 
-0.227 ±0.067 
-0.260 ± 0.087 



16 
20 
24 
32 
36 



192 
192 
192 
192 
192 



26 
24 
22 
20 
18 



0.939 
0.658 
0.595 
0.535 
0.546 



0.181 
0.539 
0.595 
0.636 
0.546 



2.641 ±0.011 
2.626 ±0.013 
2.624 ±0.019 
2.600 ±0.017 
2.597 ±0.017 



0.425 
0.422 
0.379 
0.480 
0.561 



± 0.030 
± 0.049 
± 0.074 
± 0.095 
±0.118 



0.25638 
0.25630 
0.25536 
0.25718 
0.25816 



± 0.00071 
± 0.00108 
± 0.00188 
± 0.00152 
± 0.00136 



0.2734 
0.2721 
0.2492 
0.3214 
0.4051 



±0.0132 
± 0.0235 
± 0.0333 
± 0.0743 
±0.1336 



-0.3027 ± 0, 
-0.2975 ± 0, 
-0.2222 ±0, 
-0.4843 ± 0, 
-0.9086 ± 0, 



0492 
0846 
1054 
3218 
7892 



0.566 ±0.013 
0.550 ±0.016 
0.550 ± 0.025 
0.521 ±0.018 
0.516 ±0.017 



-0.367 ± 0.034 
-0.308 ± 0.042 
-0.275 ±0.057 
-0.232 ± 0.070 
-0.266 ± 0.092 



16 
16 
16 



384 
256 
192 



1.318 
0.877 
0.914 



0.012 
0.219 
0.163 



2.612 ±0.013 
2.621 ±0.016 
2.624 ± 0.020 



0.499 
0.527 
0.516 



± 0.090 
±0.103 
±0.118 



0.25722 
0.25769 
0.25754 



±0.00115 
± 0.00126 
± 0.00149 



0.3614 
0.3854 
0.3730 



±0.1024 
±0.1199 
±0.1330 



-0.7728 ± 0, 
-0.9358 ±0, 
-0.8515 ±0, 



6469 0.6350 ± 1.0337 
8326 0.8907 ± 1.4676 
8845 0.7481 ± 1.4785 



0.533 ±0.016 
0.540 ±0.019 
0.543 ± 0.024 



-0.316 ±0.029 
-0.368 ± 0.033 
-0.372 ± 0.034 



TABLE III. List of fitting parameters corresponding to tlie global minimum Xmin of the cost function, obtained by method II for < x < 0.05. Mmin and Mmax 
represent the minimum and maximum circumferences, respectively, used in the analysis. and Q indicate the total number of input data points and goodness of fit, 
respectively. 



Mmin Mmax A'' 


xliJN Q 


f \y\ Tc/tt = ao — d 


ffli a2 as 


1 ai 


32 384 97 
36 384 86 
40 384 75 
48 384 64 
64 384 53 


1.188 0.044 
1.067 0.173 
0.631 0.978 
0.680 0.921 
0.660 0.902 


2.584 ±0.054 0.157 ±0.007 0.24629 ± 0.00084 
2.575 ±0.058 0.173 ±0.008 0.24806 ± 0.00082 
2.586 ± 0.056 0.202 ± 0.009 0.25057 ± 0.00070 
2.601 ±0.062 0.217 ±0.013 0.25150 ± 0.00084 
2.729 ± 0.060 0.260 ± 0.030 0.25360 ± 0.00126 


0.1544 ±0.0029 

0.1491 ± 0.0024 - - 
0.1437 ±0.0014 - - 
0.1428 ± 0.0010 
0.1451 ± 0.0029 


0.516 ±0.089 -0.084 ±0.192 
0.498 ± 0.089 -0.039 ± 0.214 
0.509 ±0.080 -0.063 ±0.199 
0.528 ±0.088 -0.110 ±0.217 
0.713 ±0.088 -0.529 ±0.168 


32 256 91 
36 256 80 
40 256 69 
48 256 58 
64 256 47 


1.236 0.025 
1.102 0.125 
0.595 0.986 
0.606 0.965 
0.522 0.981 


2.503 ±0.110 0.152 ±0.008 0.24572 ± 0.00103 
2.447 ±0.156 0.168 ±0.009 0.24756 ± 0.00100 
2.447 ±0.143 0.200 ±0.011 0.25037 ± 0.00085 
2.405 ±0.212 0.215 ±0.016 0.25143 ± 0.00102 
2.743 ±0.105 0.333 ±0.044 0.25606 ± 0.00112 


0.1564 ±0.0036 
0.1505 ± 0.0031 
0.1439 ±0.0018 
0.1425 ± 0.0013 
0.1565 ±0.0103 


0.391 ±0.164 0.258 ±0.580 
0.317 ±0.204 0.598 ±1.092 
0.329 ±0.171 0.584 ±0.908 
0.286 ±0.231 0.878 ±1.630 
0.705 ±0.138 -0.631 ±0.334 


32 192 85 
40 192 63 
64 192 41 


1.241 0.025 
0.583 0.983 
0.514 0.970 


2.494 ±0.130 0.149 ±0.008 0.24523 ± 0.00113 
2.392 ±0.239 0.196 ±0.012 0.25005 ± 0.00094 
2.760 ±0.131 0.342 ±0.050 0.25629 ± 0.00121 


0.1582 ± 0.0041 
0.1446 ± 0.0021 
0.1586 ± 0.0124 


0.377 ±0.193 0.305 ±0.725 
0.264 ± 0.264 0.993 ± 2.075 
0.726 ±0.171 -0.690 ±0.409 


16 384 130 
20 384 119 
24 384 108 
32 384 97 
36 384 86 


0.660 0.996 
0.689 0.985 
0.716 0.962 
0.728 0.935 
0.795 0.797 


2.655 ± 0.013 0.428 ± 0.008 0.25641 ± 0.00021 
2.644 ±0.015 0.434 ±0.012 0.25653 ± 0.00027 
2.632 ±0.019 0.412 ±0.017 0.25608 ± 0.00039 
2.629 ± 0.020 0.471 ± 0.027 0.25709 ± 0.00044 
2.620 ± 0.022 0.482 ± 0.035 0.25724 ± 0.00053 


0.2766 ± 0.0033 -0.3119 ± 0.0131 
0.2800 ± 0.0058 -0.3244 ± 0.0220 
0.2670 ± 0.0086 -0.2792 ± 0.0297 
0.3102 ±0.0199 -0.4386 ±0.0821 
0.3202 ± 0.0291 -0.4801 ± 0.1244 


0.578 ± 0.014 -0.400 ± 0.048 
0.565 ±0.016 -0.357 ±0.059 
0.553 ± 0.020 -0.294 ± 0.074 
0.546 ±0.020 -0.315 ±0.098 
0.536 ±0.023 -0.274 ±0.119 


16 256 124 
20 256 113 
24 256 102 
32 256 91 
36 256 80 


0.621 0.998 
0.631 0.996 
0.667 0.983 
0.573 0.998 
0.551 0.997 


2.655 ± 0.016 0.436 ± 0.009 0.25662 ± 0.00022 
2.640 ± 0.018 0.448 ± 0.013 0.25684 ± 0.00027 
2.623 ± 0.023 0.430 ± 0.018 0.25650 ± 0.00039 
2.612 ± 0.025 0.526 ± 0.027 0.25793 ± 0.00038 
2.596 ± 0.027 0.582 ± 0.035 0.25853 ± 0.00038 


0.2790 ± 0.0035 -0.3231 ± 0.0142 
0.2858 ± 0.0063 -0.3485 ± 0.0248 
0.2757 ± 0.0095 -0.3115 ± 0.0346 
0.3530 ± 0.0241 -0.6359 ± 0.1197 
0.4158 ± 0.0395 -0.9843 ± 0.2447 


0.577 ±0.016 -0.403 ±0.056 
0.560 ±0.019 -0.350 ±0.072 
0.542 ± 0.024 -0.268 ± 0.094 
0.527 ±0.024 -0.270 ±0.135 
0.510 ±0.025 -0.200 ±0.177 


16 192 118 
20 192 107 
24 192 96 
32 192 85 
36 192 74 


0.621 0.998 
0.634 0.994 
0.673 0.976 
0.572 0.996 
0.536 0.997 


2.658 ± 0.017 0.436 ± 0.009 0.25661 ± 0.00023 
2.642 ± 0.020 0.447 ± 0.013 0.25683 ± 0.00029 
2.623 ± 0.025 0.429 ± 0.019 0.25647 ± 0.00041 
2.610 ± 0.027 0.529 ± 0.028 0.25797 ± 0.00039 
2.590 ± 0.029 0.592 ± 0.036 0.25864 ± 0.00039 


0.2788 ± 0.0036 -0.3222 ± 0.0145 
0.2855 ± 0.0065 -0.3475 ± 0.0255 
0.2749 ± 0.0097 -0.3086 ± 0.0354 
0.3557 ±0.0251 -0.6494 ± 0.1263 
0.4273 ± 0.0421 -1.0573 ± 0.2710 


0.579 ±0.017 -0.410 ±0.058 
0.561 ±0.020 -0.355 ±0.075 
0.543 ± 0.026 -0.269 ± 0.099 
0.525 ±0.026 -0.259 ±0.144 
0.505 ±0.027 -0.169 ±0.193 


16 384 130 
16 256 124 
16 192 118 


0.657 0.995 
0.583 0.999 
0.584 0.999 


2.651 ± 0.013 0.461 ± 0.045 0.25683 ± 0.00059 

2.644 ±0.014 0.529 ±0.033 0.25776 ± 0.00040 

2.645 ± 0.015 0.530 ± 0.034 0.25777 ± 0.00042 


0.3122 ± 0.0487 -0.4850 ± 0.2548 0.1989 ± 0.3207 
0.3829 ± 0.0382 -0.9177 ± 0.2648 0.8419 ± 0.4628 
0.3832 ± 0.0395 -0.9203 ± 0.2749 0.8461 ± 0.4808 


0.570 ±0.015 -0.410 ±0.052 
0.559 ±0.014 -0.438 ±0.062 
0.560 ± 0.014 -0.443 ± 0.064 



TABLE IV. List of fitting parameters corresponding to the global minimum Xmin of the cost function, obtained by method II for < a; < 0.07. Mmin and Mmax 
represent the minimum and maximum circumferences, respectively, used in the analysis. N and Q indicate the total number of input data points and goodness of fit, 
respectively. 



Mmin Mmax N 


xlnJN Q 


f \y\ Fc/tt = ao — d 


ai a2 as 


1 «! 


32 384 131 
36 384 116 
40 384 101 
48 384 86 
64 384 71 


1.225 0.018 
1.214 0.025 
0.904 0.589 
0.936 0.463 
1.066 0.172 


2.588 ±0.025 0.146 ±0.006 0.24477 ± 0.00085 
2.595 ± 0.026 0.160 ± 0.007 0.24655 ± 0.00083 
2.606 ±0.025 0.188 ±0.008 0.24934 ± 0.00071 
2.593 ±0.030 0.203 ±0.011 0.25045 ± 0.00083 
2.617 ±0.036 0.218 ±0.026 0.25143 ± 0.00154 


0.1601 ± 0.0031 
0.1543 ±0.0027 

0.1471 ± 0.0017 - - 
0.1451 ± 0.0014 - - 
0.1446 ± 0.0009 


0.535 ±0.043 -0.134 ±0.086 
0.542 ±0.043 -0.150 ±0.085 
0.553 ±0.038 -0.180 ±0.080 
0.529 ±0.044 -0.128 ±0.103 
0.563 ±0.052 -0.212 ±0.121 


32 256 123 
36 256 108 
40 256 93 
48 256 78 
64 256 63 


1.227 0.019 
1.211 0.029 
0.844 0.731 
0.840 0.692 
0.866 0.567 


2.600 ± 0.034 0.147 ± 0.007 0.24495 ± 0.00096 
2.611 ± 0.036 0.164 ± 0.008 0.24696 ± 0.00093 
2.627 ±0.036 0.197 ±0.010 0.25005 ± 0.00077 
2.602 ±0.048 0.216 ±0.014 0.25138 ± 0.00090 
2.691 ±0.051 0.307 ±0.038 0.25529 ± 0.00115 


0.1595 ± 0.0035 
0.1529 ±0.0030 
0.1454 ±0.0017 
0.1436 ±0.0011 
0.1520 ±0.0072 


0.555 ±0.059 -0.171 ±0.108 
0.567 ±0.058 -0.197 ±0.108 
0.581 ±0.053 -0.237 ±0.103 
0.539 ±0.067 -0.154 ±0.156 
0.646 ±0.065 -0.469 ±0.171 


32 192 115 
36 192 100 
40 192 85 
48 192 70 
64 192 55 


1.217 0.024 
1.211 0.031 
0.827 0.746 
0.830 0.684 
0.841 0.585 


2.606 ± 0.037 0.145 ± 0.007 0.24461 ± 0.00104 
2.618 ±0.039 0.162 ±0.008 0.24673 ± 0.00101 
2.636 ±0.039 0.196 ±0.010 0.24999 ± 0.00083 

2.607 ±0.054 0.215 ±0.015 0.25133 ± 0.00098 
2.731 ±0.060 0.327 ±0.043 0.25585 ± 0.00116 


0.1608 ± 0.0039 
0.1537 ± 0.0033 
0.1456 ±0.0018 

0.1437 ±0.0012 - - 
0.1559 ± 0.0097 


0.567 ±0.064 -0.191 ±0.113 
0.580 ±0.064 -0.220 ±0.113 
0.595 ±0.058 -0.262 ±0.109 
0.547 ±0.077 -0.172 ±0.173 
0.694 ± 0.076 -0.601 ± 0.200 


16 384 176 
20 384 161 
24 384 146 
32 384 131 
36 384 116 


0.820 0.915 
0.810 0.916 
0.830 0.860 
0.873 0.721 
0.957 0.429 


2.624 ±0.006 0.425 ±0.007 0.25623 ± 0.00019 
2.612 ± 0.007 0.425 ± 0.010 0.25623 ± 0.00025 
2.604 ±0.009 0.400 ±0.015 0.25567 ± 0.00037 
2.595 ± 0.010 0.438 ± 0.024 0.25640 ± 0.00046 
2.597 ±0.010 0.464 ±0.031 0.25682 ± 0.00051 


0.2791 ± 0.0028 -0.3149 ± 0.0111 
0.2798 ± 0.0048 -0.3173 ± 0.0183 
0.2655 ± 0.0071 -0.2678 ± 0.0243 
0.2921 ± 0.0163 -0.3602 ±0.0611 
0.3127 ± 0.0240 -0.4402 ± 0.0983 


0.549 ± 0.006 -0.317 ± 0.023 
0.536 ± 0.007 -0.266 ± 0.029 
0.528 ±0.009 -0.219 ±0.036 
0.518 ±0.010 -0.195 ±0.047 
0.519 ±0.011 -0.216 ±0.056 


16 256 168 
20 256 153 
24 256 138 
32 256 123 
36 256 108 


0.706 0.995 
0.691 0.995 
0.727 0.979 
0.702 0.982 
0.704 0.970 


2.633 ± 0.007 0.436 ± 0.008 0.25655 ± 0.00019 
2.619 ±0.008 0.443 ±0.011 0.25666 ± 0.00025 
2.610 ±0.011 0.423 ±0.016 0.25626 ± 0.00035 

2.598 ± 0.012 0.497 ± 0.025 0.25745 ± 0.00038 

2.599 ± 0.012 0.566 ± 0.030 0.25829 ± 0.00036 


0.2819 ± 0.0030 -0.3297 ±0.0122 
0.2862 ± 0.0053 -0.3454 ± 0.0209 
0.2749 ± 0.0080 -0.3040 ± 0.0288 
0.3316 ±0.0199 -0.5285 ± 0.0900 
0.4021 ± 0.0330 -0.8916 ±0.1936 


0.557 ±0.007 -0.349 ±0.026 
0.541 ± 0.008 -0.296 ± 0.034 
0.534 ±0.011 -0.250 ±0.043 
0.518 ±0.011 -0.229 ±0.062 
0.518 ±0.012 -0.282 ±0.080 


16 192 160 
20 192 145 
24 192 130 
32 192 115 
36 192 100 


0.690 0.996 
0.677 0.996 
0.714 0.980 
0.686 0.984 
0.677 0.978 


2.635 ± 0.007 0.437 ± 0.008 0.25656 ± 0.00020 
2.621 ±0.009 0.443 ±0.011 0.25667 ± 0.00025 
2.612 ±0.011 0.422 ±0.016 0.25625 ± 0.00037 

2.598 ± 0.013 0.500 ± 0.026 0.25751 ± 0.00040 

2.599 ± 0.013 0.577 ± 0.031 0.25843 ± 0.00036 


0.2819 ± 0.0030 -0.3299 ±0.0125 
0.2863 ± 0.0054 -0.3460 ± 0.0215 
0.2746 ± 0.0082 -0.3031 ± 0.0297 
0.3344 ± 0.0208 -0.5415 ± 0.0954 
0.4139 ± 0.0352 -0.9636 ± 0.2158 


0.559 ± 0.007 -0.356 ± 0.027 
0.543 ± 0.009 -0.302 ± 0.035 
0.535 ±0.011 -0.255 ±0.045 
0.518 ±0.012 -0.231 ±0.065 
0.517 ±0.012 -0.287 ±0.086 


16 384 176 
16 256 168 
16 192 160 


0.817 0.911 
0.680 0.998 
0.664 0.998 


2.620 ± 0.007 0.456 ± 0.039 0.25663 ± 0.00053 
2.622 ± 0.007 0.516 ± 0.030 0.25753 ± 0.00038 
2.624 ±0.007 0.517 ±0.032 0.25756 ± 0.00040 


0.3122 ± 0.0419 -0.4751 ± 0.2170 0.1817 ± 0.2679 
0.3704 ± 0.0346 -0.8216 ± 0.2280 0.6671 ± 0.3728 
0.3710 ± 0.0359 -0.8267 ± 0.2372 0.6746 ± 0.3891 


0.543 ± 0.009 -0.322 ± 0.025 
0.541 ± 0.007 -0.364 ± 0.030 
0.543 ±0.008 -0.371 ±0.031 
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TABLE V. In this table we compare our results for the exponents v, \y\, Vc, and ao — 2 with those of other authors. In the 
second column the abbreviations "CC model" , "LLL" , and "TBM" stand for the Chalker-Coddington network model, the lowest 
Landau level basis approximation, and the tight-binding model in a perpendicular magnetic field, respectively. Also, systems 
with different dimensionality are denotes as "2D" (two dimensions), "QID" (quasi-one dimension), and "ID" (one dimension). 
"PRC" and "RBC" stand for the periodic and reflecting boundary conditions, respectively, in the transverse direction in QID. 



model 



present work 



CC model with PBC in QID 



2.62 ±0.06 



>0.4 



0.257 ±0.002 



Huckestein et al. 
Mieck g 
Huckestein Q 
Huo et al.^ 
Huckestein ^ 
Lee et alM 
Cain et al.\7\ 
Sandler et al. 
Slevin et al. \ 
Obuse et al. 
Amado et al. ]f 1 | 
Dahlhaus et al. |^ 
Fulga et al. |I3] 
Slevin et al. [|l4| 
Janssen 
Evers et al. |I( 
Obuse et al. |f 7| ] 
Evers et al. ill 



LLL in QID 
LLL in QID 
LLL in QID 
TBM in 2D 
LLL in QID 
CC model in QID 
CC model in 2D 
LLL in 2D 

CC model with PBC in QID 

CC model with RBC in QID 

CC model with PBC in QID 

kicked rotator in ID 

CC model in Corbino disc 

CC model with PBC in QID 

CC model in 2D 

CC model in 2D 

CC model in 2D 

CC model in 2D 



2.34 ±0.04 

2.3 ±0.08 
2.33 ±0.05 

2.4 ±0.1 

2.33 ±0.03 
2.37 ±0.02 
2.33 ±0.09 
2.593 ± 0.006 

2.55 ±0.01 
2.616 ±0.014 
2.576 ± 0.03 

2.56 ±0.03 
2.607 ± 0.009 



0.38 ±0.04 



0.17 
1.29 

(logarithmic correction) 
no corrections 



0.279 ±0.05 



0.248 ± 0.004 



0.223 ±0.05 



0.27 ±0.02 
0.262 ± 0.003 
0.2617 ±0.0006 
0.2596 ± 0.0004 



